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Abstract. In the present paper we propose a new stochastic diffusion process with 
drift proportional to the Weibull density function defined as 



t > 0, 


with parameters 7 > 0 and ct > 0, where B is & standard Brownian motion and t = e 
is a time near to zero. First we interested to probabilistic solution of this process 
as the explicit expression of this process. By using the maximum likelihood method 
and by considering a discrete sampling of the sample of the new process we estimate 
the parameters 7 and <j. 
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1 Introduction 

In the present paper we propose a new stochastic Weibull process X = {Xt, t > 0} 
given by the following linear stochastic differential equation 



t > 0 


where /i(f, W;7) = — CW and g{Xt]a) = aXt. We denote by i? a 


standard Brownian motion and 7 and a are an unknown parameters. An interesting 
problem is to estimate the parameters 7 and a are time independent reel parameters 
to be estimate when one observes the whole trajectory of X. 

The estimation for diffusion processes by discrete observation has been studied by 
several authors ( see for example [S], [Z], PP, [2], [S])and its references. Prakasa Rao 
[7] treats this problem and shows that the least square estimator is asymptotically 
normal and efficient under the assumption h\fN —)■ 0, the condition for ’’rapidly 
increasing experimental design” [7]. 

The organization of our paper is as follows. Section 2 contains the presentation 
of the basic tools that we will need throughout the paper: basic properties of of 
standard Brownian motion and It’s formula . The aim of section 3 is twofold. 
Firstly, we prove the close formula of SDE under the conditions (1) — (3). Secondly, 
we invesigate the mean of Xt an explicit solution of SDE and we prove that the drift 
is proportionnel to Weibul density. The section 4 is devoted to estimate a parameters 
by using Maximum likelihood. In the last section we present a numerical test. 
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Given the general one-dimensional time-homogeneous SDE 

Xe = x > 0 ; dXt = fi{Xt]'y)dt + g{Xt]a)dBt, t > 0. ( 1 . 1 ) 

where fi{t, Xt] 7 ) = y(l — — PXt and g{Xt] a) = aXt. The following condi¬ 

tions are assumed in this article : 

(1) There is a constant L > 0 such that 

|/i(t,x; 7 )| + \g{x]a)\ < L (1 + |x|) 

(2) There is a constant L > 0 such that 

|/i(t,x; 7 ) -/i(t,|/; 7 )| -F \g{x]a) - g{y]a)\ < L\x - y\ 

(3) For each g > 0, sup^E(|Xt|'^) < 00 . 

2 Preliminaries 

In this section we start by recalling the dehnition of Brownian motion, which is a 
funda- mental example of a stochastic process. The underlying probability space 
(G, P) of Brownian motion can be constructed on the space G = Co(M’'") of con¬ 
tinuous real-valued functions on M"*" started at 0. For more complete presentation 
on the subject, see a, 0 . 

Definition 2.1. The standard Brownian motion is a stochastic process {Bt, t > 0) 
such that 

(a) Bq almost surely; 

(b) With probability one t —)■ Bt is continuous. 

(c) For any finite sequence of times t^ < ti < ■ ■ ■ < t^ the increments 

Bit — BtQ, Bt^ — Btj^, ■ ■ ■ Bt^ — Bt^_^ 

are independent. 

(d) for any given times 0 < s < f, Bt — Bg has the Gaussien distribution N, t—s) 
with mean zero and variance t — s. 

We refer to Theorem 10.8 0 / 0 / and to Chapter of 0/ for the proof of the existence of 
Brownian motion as a stochastic process {Bt, t > 0) satisfying the above properties 
(a) — (d). In the sequel the filtration {iFt)t>o will be generated by the Brownian paths 
up to time t, in other words we write 

Tt = (y{Bs : 0 < s < f), f > 0. 

we give a basic properties of Brownian motion and extentions ([?]): 

• The crucial fact about Brownian motion, which we need is {dBfi = dt] 

• For every 0 < s < f, Bt — Bg is independent of {B^, u < s} and has a 
Af{0,t-s); 
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• Brownian motion [Bt) is a process Markov property; 

• {—Bt)t>o is a Brownian motion; 

• Xt = Xq + /it + aBt is a Brownian motion with drift with mean eqnal to 
Xq + fit] 

• Xt = Xoexp(/it + aBt) is a Geometric Brownian motion with mean eqnal to 
Xq exp(/it + a"^/2). 

We introduce the following two version of It’s formula 
Theorem 2.2. ( ltd’s formula v.l) Let f G Then for a < t, 

= r{B.) * + i r !"iB.) ds (2.1) 

Theorem 2.3. ( ltd’s formula v.2) Let = f{t,x) be a continuous in [a, 6] x M with 
ft, fx, and fxx continuos (a, b) x M. Then for a < t < b, 

f{t,Bt)-f{a,B^) = Jjx{s,B,)dB, + J^(^ft{s,B,) + ^fxx{s,B,)^ ds. (2.2) 


3 Closed formula and mean of X 


t 


By curiosity, we focus on the explicit formula for the SDE to hnd the mean and 
variance of Xt. By using the It’s formula to h) = In(Xt) we obtain 


dYt = ( |(1 - - y ) + (^dBt, W = Inx. 


a 


(3.1) 


By integrating between e and t it follows 

t 


Tt = In x + 7 In 


(r+' - £'*+') - y); - e) + <t(B, - B, 


(3.2) 


Then the explicit solution of SDE is given by 


Xt = x^^ exp - e) + (T(5t - , Vf > 0. (3.3) 

We intersted now to mean of Xt. By using the conditional expectation or the 
Geomtric Brownian we prove that the trend of the process Xt is given by 


E(W) = X 


exp (e''+‘ - r+‘), 


(3.4) 


then it follows that the trend of Xt is proportional to Weibul density. 
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4 Maximum likelihood estimators for the param¬ 
eters of SDE 


A formal statement of the parameter estimation problem to be addressed is as 
follows. Given the general one-dimensional time-homogeneous SDE 


ATe = X > 0; dXt = + g{Xt;a)dBt, t > 0. (4.1) 

where n{t,Xt;'y) = 7(1 - t'^^^)Xt - PXt and g{Xt;a) = aXt. 

The task is to estimate the parameters 6 = ( 7 , a^) of this SDE from a sample 
of TV -|- 1 observations Xi, X 2 , ■ ■ ■ , X^ oi the stochastic process at known times 
ti, ■ ■ ■ ,tN ■ In the statement of equation fl4.ip . dB is the differential of the Brownian 
motion and the instantaneous drift fi{t,x]'y) and instantaneous diffusion g{x]a). 
Assuming that = x) = p ^ 0, for the sake of simplicity, let us assume that 

p = 1. 

The ML estimate of 6 is generated by minimising the negative log-likelihood 
function of the observed sample, namely 

N 

ll{x,,X 2 ,--- ,x^-e) = iogMx,\e) -J2^ogf{Xk+i\Xk,e) (4.2) 

k=2 


with respect to the parameters 9 = ( 7 , a^). In this expression, /i(JVi| 6 *) is the density 
of the initial state and f{Xk+i\Xk; 6) = /((X^+i, tfe+i)|(JVfc, tk)] 6) is the value of the 
transitional PDF at (JVfc+i,tfc+i) for a process starting at {Xk,tk) and evolving to 
{Xk+i,tk+i) in accordance with equation fl4.ip . Note that the Markovian property 
of equation fl4.ip ensures that the transitional density of Xk+i at time tk+i depends 
on Xk alone. 

ML estimation relies on the fact that the transitional PDF, f{x, t), is the solution 
of the Fokker-Planck equation 


dt 


^ ( l g(^(x;a)/) 
\2 dx 




(4.3) 


satisfying a suitable initial condition and boundary conditions. Suppose, further¬ 
more, that the state space of the problem is [a, b] and the process starts a.t x = Xk 
at time tk- In the absence of measurement error, the initial condition is 


f{x,tk) = S{x - Xk) 


(4.4) 


where 6 is the Dirac delta function, and the boundary conditions required to conserve 
unit density within this interval are 


lim 

X — 


l d{g{x; a)f) 
2 dx 


= 0 , 


lim , 

X — yb~ V 2 


ld{g{x;a)f) 


dx 


- ) = 0 . 


(4.5) 


By using Ito formula to L) = In(Xj) we have the linear diffusion 
dYt = (|(1 - - y ) dt + (jdBt, 

By integrating between tk and tk+i, the exact discret model correspond to fl4.14p is 
given by 
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From the above, the variance and esperance of ln(Xfc+i) is given by 


E(ln(Xt+,)l Vt) = In(Vt) + 71.1 - ((2:1 - *2+') - y (4+i - h), 

V (ln(Xfc+i) \Xk) = cr^(tfc+i — tfc), 


then the transitional probability density function (PDF) for SDE has the following 
closed form expression: 

Xt+i\Xt ~ Af fln(Jfi) + 7 In - (* 22 i - * 2 ^') - “ 4 ),ct^( 4 +i - 4 )f 

(4.7) 

Now, the classical approach to the ML method requires the computation of the 
hrst-order partial derivatives of the log-likelihood function with respect to each of 
its parameters, equating them equal to zero and then solving the resulting system 
of equations. So, the hrst-order partial derivatives are obtained as follows: 


dLL 

d'y 

dLL 


5^E(''' + 2-4»(W.4))(ln(^)-(7 + l)(i2-(2-i)) = 0 . (4.8) 

- 11 ^ 1 ^ 

“ 2^2 + 3^ + 274fc(Xfc, + ‘^^k{Xk, tk)) = 0 

k=2 k=2 

(4.9) 


with Ak{Xk,tk) = In(^) - tlt\. 

From the equation fl4.9l) we obtain 


N 


i(N - l)a^ - (N - l)a* = iY, 4) 


(4.10) 


k=2 


By using the positivity of we have the following expression of estimator 


N 


1/2 


F-ivVTE4‘2(W.4)) - 2 . 


(4.11) 


k=2 


Consequently, we prove the non-linear expression of estimator 7 by replacing d^ in 

(gSl) : 

f^(ff= + 2i»(At,4 ))(^ln(^)-(7+l)(i2-i2 -i)) =0 (4.12) 

where Ak{Xk,tk) = In(^) - 7 ln(^) + 

It is assumed that the observation from the realization consists of ti = 
e, tk = kh with d > 0, k = 2,3, ■ ■ ■ ,N. We dehne the new estimator of : 


d^ = 4 


X- 1 


^ / V u \ ■ 

^ ^ d^+^(F+^ -{k- 1)^+^) j 


1/2 


- 2 . 


(4.13) 


Assuming that 7 = 7. The following proposition is more or less well khnown. 
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Proposition 4.1. If condition (1) — (3) hold, then 

E(ff^ -a^)<C + ft) j . 


Proof. By using integrating the following equation between t^-i = {k — l)h and 
tk = kh\ 

2 \ 

— P—Yjdt + adBt, (4-14) 

By summing between k = 2 and N and using = h. Then we prove 

a result. □ 


dYt= (^2(1-r+i) 


5 Numerical test 

We present numerical results in the software R for the following SDE: 

X, = 10, dXt = (|(1 - Xtdt + aXtdBt, t > 0, 

with: and we generate sampled data with 7 = 1, a = 0.2 and time step A = 10“^ 
as following: 

• mu <- expression( ((l./t)- 2 *t)*x) 

• g <- expression( 0 . 2 *x) 

• Simulation <- SNSDE(drift=mu,diffus=g,N=1000,Dt=0.001,x0=10) 

• MyData <- SimulationX 

• mux <- expression( gama * x./t — (gama + 1 ) * fiaama+i) ^ ^ ^ 

• gx <- expression ( sigma * x ) 

• Model <- FDE(data=MyData,drift=mux,diffus=gx,start = hst(gama=l,sigma=0.2),”euler”) 

• summary (Mo del) 

Result : 

Pseudo maximum likelihood estimation 

Method: Euler 
Call: 

FSDE(data = MyData, drift = mux, diffus = gx, start = list(gamma=l, sigma=l),”euler”) 

The following table prove the estimate coefficients and standard error with —2 log L : 

8794.337 by unsig our methods. 



Estimate coefficient 

Standard error 

7 

a 

0.9856006 

0.2310434 

0.006651509 

0.005168557 


Figure 1: Estimate coefficients and Standard Error 


By using the command in the R software ” conhnt(Model,level= 0.9)”, we obtain 
the following table. It shows the conhdence interval for the estimated variables 7 
and a. 
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5% 

15% 

75% 

95% 

7 

a 

0.9711487 

0.2225419 

1.0000124 

0.2238546 

1.0003527 

0.2310215 

1.0012685 

0.2395449 


Figure 2: Confidence interval 
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